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We analyze the Mott-insulator phases of dipolar bosonic gases placed in neighboring but uncon- 
nected ID traps. Whereas for short-range interactions the ID systems are independent, the non-local 
dipole-dipole interaction induces a direct Mott-insulator to pair-superfluid transition which signif- 
icantly modifies the boundaries of the lowest Mott-insulator phases. The lowest boundary of the 
lowest Mott regions becomes progressively constant as a function of the hopping rate, eventually 
inverting its slope, leading to a re-entrant configuration which is retained in 2D. We discuss the 
consequences of this effect on the spatial Mott-insulator plateaux in experiments with additional 
harmonic confinement, showing that anti-intuitively the plateaux may become wider for increasing 
hopping. Our results are also applicable to non-dipolar boson-boson mixtures. 



Strongly-correlated atomic gases have recently at- 
tracted a rapidly-growing attention, mainly motivated by 
impressive developments on the manipulation of atoms 
in optical lattices. When loaded in these lattices, ultra 
cold atoms experience a periodic potential that resembles 
that of electrons in solids, opening fascinating links be- 
tween the physics of cold atoms and solid-state physics. 
In particular, cold bosons restricted to the lowest lattice 
band can be described by the Bose-Hubbard model 
which presents two different phases at zero temperature 
Q, namely a superfluid (SF) phase, and a gapped incom- 
pressible insulator phase known as Mott-insulator (MI), 
characterized by a commensurate occupation per lattice 
site. The SF to MI transition in cold bosons in optical 
lattices was recently observed in a remarkable experiment 
Q , in which the gapped nature of the MI excitation spec- 
trum was clearly demonstrated. The realization of the MI 
was indeed possible due to an additional harmonic poten- 
tial overimposed to the lattice that guaranteed locally the 
(otherwise practically impossible) commesurability con- 
dition necessary for the MI. Such inhomogeneous poten- 
tial leads to the formation of spatial MI and SF shells 
EL 01 j which have been observed very recently 0, @| ■ 

Among the various research lines related with optical 
lattices, the physics of mixtures has attracted a consider- 
able attention. Bose-Fermi mixtures may lead to a wealth 
of phases of fermion composites 0, and may allow for 
the generation and engineering of disorder [1, 0, [l^, 11 1 . 
Bose-Bosc mixtures have also attracted a major interest 
l|, Q m- I- Particular, it has been shown 
1 31 ] that the interspecies interaction may lead to the for- 
mation of a pair superfluid (PSF), i.e. a superfluid of 
boson-boson (or hole-hole) composites, which occurs in 
addition to MI phases for both components, as well as 
uncorrelated superfluid phases in each one of the compo- 
nents (2SF). 

Dipolar gases also attract currently a major interest, 
motivated by recent experiments on atoms with large 

polar molecules [18| and Ry- 
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nificant or even dominant when compared to the short- 
range isotropic interactions [jo| . The DDI can play an 
important role in the physics of lattice bosons, leading to 
additional phases, as checker-board or supersolid, which 
may be easil y c ontrollable by manipulating the atomic 



confinement 2l|. In addition, contrary to short-range 



interacting gases for which disconnected sites (i.e. with- 
out hopping between them) are fully independent, the 
long-range character of the DDI induces a coupling even 
for unconnected sites. The latter leads to fundamentally 
new physics, as e.g. a condensate of filaments novel 
quantum pha se transitions in bilayer systems of polar 
molecules 23], or inelastic interlayer scattering for dipo- 
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magnetic moment 

dberg atoms [lij]. In these gases the long-range and 
anisotropic dipole-dipole interactions (DDI) become sig- 



lar solitons 

In this Letter, we show that the non-local DDI in- 
duces a direct MI-to-PSF transition for neighboring un- 
connected ID dipolar gases, which significantly modifies 
the boundaries of the lowest ID MI phases. The same 
effect is also expected under appropriate conditions for 
non-dipolar boson-boson mixtures. Remarkably, the low- 
est boundary of the first MI lobes eventually inverts its 
slope as a function of the hopping rate, leading to a re- 
entrant scenario, which is maintained in 2D. We show 
that this effect has direct consequences on the spatial ex- 
tension of the MI plateaux for the case of an overimposed 
harmonic confinement. Anti-intuitively, we show that the 
plateaux extension may become constant or even wider 
for increasing hopping. 
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Figure 1: Scheme of the system considered in this Letter. 

In the following, we consider dipolar bosons placed 
at two neighboring, but disconnected, ID traps (wires), 
which can be created using micro-magnetic confinement 
j2f| or sufficiently strong 2D optical lattices. In the latter 
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case, the required two-site configuration may be gener- 
ated by super-lattice techniques or by selectively empty- 
ing ID sites neighboring the desired pair. Along the ID 
systems we assume an additional lattice equal for both 
ID traps, which leads to the ladder configuration shown 
in Fig. [T] In the following we borrow from the literature 
on spin ladders the terms leg (each wire) and rung (pair 
of neighboring sites belonging to different legs). In this 
Letter, we are mostly concerned about interlayer effects, 
and hence we consider a configuration for which only the 
(attractive) DDI between sites at the same rung plays a 
significant role. This is the case, if the dipoles are ori- 
ented forming an angle <p with the axis of the wires, such 
that cos 2 (p = 1/3. In that case, the DDI between neigh- 
bors at the same leg vanishes, whereas the DDI between 
sites in the same rung is attractive. There is in princi- 
ple an additional non zero diagonal DDI between sites 
in neighboring rungs belonging to different legs. These 
terms can be made negligible by considering the spacing 
between rungs, 7 > 1 times larger than the separation 
between the two legs. In that case the spurious diagonal 
interaction is a factor (1 + 2y/2j)/(l +7 2 ) 5/2 (~ 0.03 for 
7 = 3) smaller than that between sites in the same rung. 
Of course, for other dipole and lattice configurations, the 
DDI between sites belonging to the same leg cannot be 
neglected, and interesting physics can be expected [26| 
and will be studied elsewhere. 

Under the previous conditions the system is described 
by a Bose-Hubbard Hamiltonian (BHH) of the form 

a=l,2 <i,j> a=l,2 



where b\ a \ b[ a ^ , and nl"' are, respectively, the annihila- 
tion, creation, and number operators for the site i at the 
leg a. J describes the hopping between neighboring sites 
i and j in each leg, Uq the on-site interactions (a combi- 
nation of short-range and dipolar contributions [2l|), and 
we consider the same chemical potential [i in both legs. 
Atoms in sites at the same rung interact attractively by 
the DDI, which is characterized by a coupling — |J7'|. 

In the following we analyze the effects of the coupling 
U 1 in the physics of the MI phases for the ID wires. Note 
that the Hamiltonian ([1]) is formally equivalent to the 
case of two bosonic species in a ID array, with equal 
chemical potential fi for both, equal hopping J, equal 
on-site intra-species interactions Uq, and an interspecies 
interaction — \U'\. Hence, our results can be equally ap- 
plied to boson-boson mixtures under these constraints. 
Indeed, as we show below, the PSF phase introduced in 
the context of Boson-Boson mixtures 131 is crucial for 



the understanding of the physics discussed below. 

In our analysis of the ground states of H , we have em- 
ployed Matrix-Product-State (MPS) techniques, follow- 
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Figure 2: MI lobes as a function of J/Uo and [i/Uo, for 
\U'\/Uo = (a), 0.25 (b), 0.5 (c) and 0.75. The calcula- 
tions where made for L — 12 and D = 5. The error bars 
indicate the change in the order parameters, (l)$psp and 
(•)(&), from < 10 -4 to > 10~ 2 . 



ing closely the method of Ref. (27j. The MPS represent 
an optimal Ansatz [28| for problems as that of this Let- 
ter. Adapted to the two-leg problem, with L sites per 
leg, the MPS Ansatz for the many-body wavefunction is: 



I* >= 



E 

(o)_ 0} 



A, 



(1) (2)i 
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l \{n^}), 



where we consider a maximal number of atoms n r 



(2) 



per 



site, and A 



is a D x D matrix, associated to the 



case of rij and rij atoms at the site j of both legs. In 
the following we focus on the regime of low average occu- 
pation per site around the first MI lobe, and hence in our 
calculations it proves enough n max — 2. In addition, we 
have checked in our calculations that relatively low ma- 
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trix dimensions D = 6 describe properly the problem un- 
der consideration. The MPS Ansatz enormously simplify 
the original problem (which scales exponentially with L), 
since it has a complexity given by {rij rinT + 1)D 2 L. Using 
a similar approach as that of Ref. [27j we developed a nu- 
merical algorithm that allows us to recursively adapt the 
matrices until reaching the ground state. This method re- 
sembles in many ways that of finite-size Density-Matrix- 
Renormalization-Group techniques [29| . 
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Figure 3: (M)ippsF and (•){&) as a function of J/Uo for n — 
-0.15 and \U'\/U = 0.75 with L — 24 and D = 6. 

Fig. [2] shows the results of our simulations for the sur- 
roundings of the lowest MI lobe (with (n ? 1,2 ) = 1) for 
\U'\/U o = (a), 1/4 (b), 1/2 (c) and 3/4 (d). Note that 
in order to avoid collapse in a single site, \U'\ < Uq. For 
the case of U' = 0, the usual Mott-lobes are recovered 
[30j | . However, the dependence of the lobe boundaries 
in the \i-J phase space changes significantly when \U'\ 
grows. Note, in particular, that the lowest boundary be- 
comes progressively flatter when \U'\ approaches Uo/2. 
Indeed our analytical calculations (see below) show that 
for \U'\ > Uo/2 the slope of the lowest boundary of the 
MI lobe inverts its sign. This behavior is however diffi- 
cult to observe in our numerical calculations due to the 
very narrow region between the Mi-lobe and the region 
of zero occupation. In the following we discuss in more 
detail the physics behind the distortion of the MI lobes, 
and the implications of this distortion on the spatial ex- 
tension of the MI lobes in axially trapped gases. 

The boundaries of the MI lobes are provided by the 
energy gap between the MI state and the lowest excited 
state conserving the particle number. In usual (single- 
component) BHH this lowest excitation is provided 
by particle-hole excitations. The MI boundaries can then 
be calculated by a strong-coupling expansion (SCE) [3l1 |. 
estimating the energy of a state with an extra particle 
and a state with an extra hole. This is indeed the case 
of U' = 0, where the lowest excitations are given by un- 
corrected particle- hole excitations in both wires. The 
situation changes for \U'\ > 0, since for sufficiently low 
tunneling, there is a direct transition between MI and 
PSF phases, i.e. supcrfluid phases of composites (or com- 
posite holes) [l]| . In that case the first excitation of the 



MI lobe is given by the correlated creation of pairs of 
particles (or holes) at opposite sites of the two wires, ex- 
plaining the qualitative change in the shape of the lobe 
boundaries. In particular, a second-order SCE in J/\U'\ 
[12], provides the following dependence for sufficiently 
low tunneling for the lowest boundary of the MI lobe 
with uq particles per site: 
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From ([3]) it becomes clear that for any U' > the gap 
boundaries are quadratic (and not linear) in J for suf- 
ficiently low J. Interestingly, the lowest boundary of 
the first MI region (no = 1) inverts its slope at J = 
for \U\' > Uo/2, in agreement with our numerical re- 
sults. One may also observe that an inversion of the 
slope of the lowest boundary is expected also for no = 2 
at \U'\/Uq ~ 0.85, but it is not expected for no > 2. 

In our numerical simulations, we revealed the presence 
of the pairing-phases by monitoring ^psf = K^ 1 **^ 2 -*) — 
(SW)(S( 2 )}| (33[. A typical dependence of * PSF and 
(b) = (M 1 - 2 )) in our simulations is depicted in Fig. [3] 
for a fixed chemical potential. The MI region is char- 
acterized as that in which both fypsF = (b) = 0. We 
denote the PSF region as that in which ^>psf ^ but 
(b) = 0. Finally the 2SF region is that in which (b) =/= 0. 
Note that there is a finite coexistence region in which 
both ^>psf 7^ and (b) ^ 0. Repeating the calculations 
for different chemical potentials we obtain the results de- 
picted in Figs. [2 Note that, as we mentioned above, a 
direct MI-PSF transition can be observed in Figs. [5] at 
low J/Uq, which results in a clear distortion of the MI 
boundaries when compared to the U' = case. 

The qualitative change in the shape of the MI lobe 
has important consequences on the spatial extension of 
the MI and SF regions in the presence of an overimposed 
harmonic confinement. In order to analyze this point, 
we consider a harmonic trap along the wires, such that a 
term Q J2i a i 2 nf^ is added to the BHH. This term in- 
duces a local chemical potential /i; = /io — £li 2 , where uq 
is the local chemical potential at the trap center. Hence, 
Hi scans values (i < /iq. If for a given tunneling rate, 
the system with is inside the first MI lobe, it is hence 
expected the apperance of a MI shell at the trap cen- 
ter, characterized by a plateau in the average popula- 
tion per site ((n) = 1), surrounded by a second SF shell 
(with (n) < 1) (H, For U' = 0, for a fixed chemical 
potential, it is intuitively expected that the MI plateau 
shrinks when J/Uq increases, until eventually dissappear. 
Indeed, this is the case, since the lowest boundary of the 
MI lobe increases with J, hence decreasing the spatial 
MI region (Fig. IUa)). However, when \U'\ grows, the 
change in the slope of the lowest boundary of the first 



4 



MI lobe leads to a significant modification of the spatial 
extension of the MI plateau. In particular, as shown in 
Figs. |3J the basically J-independent lowest MI bound- 
ary for \U'\ = Uq/2 leads to a J-independent MI plateau 
(Fig. (He)). Moreover, for \U'\ > U /2, the re-entrant 
character of the MI lobe leads to the anti-intuitive ob- 
servation, that for enhanced mobility (J larger) the MI 
plateaux become even broader (Fig. HJd)). 
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Figure 4: Spatial distribution of (n) for U' /Uo = (a), 
0.25 (b), 0.5 (c) and 0.75 (d), for, respectively 10 4 Q/U = 
12.3,9.26,9.26,4.32, with fi /U = 0.2,0,-0.1,-0.36, for 
J/Uo = 0(B), 0.05(«), 0.1(A). 

Although a detailed 2D analysis is beyond the scope 
of this Letter, and will be the subject of further inves- 
tigations, we stress here that the SCE for 2D lattices 
at unconnected layers (or equivalcntly to two-component 
bosonic gases in 2D lattices) shows that the lowest 
boundary of the first MI lobe follows at low J the re- 
lation © but substituting 2{J/U ) 2 by z(J/U ) 2 , where 
z is the coordination number. Hence, the change in the 
sign of the slope occurs exactly as for ID, and thus a 
re-entrant scenario is also expected in 2D. 

In this Letter we have analyzed the physics of dipolar 
gases in unconnected neighboring ID systems. Whereas 
without dipolar interactions the ID systems are inde- 
pendent, the nonlocal dipolc- induced intcrlaycr interac- 
tion leads to a direct MI to PSF transition, significantly 
distorting the Mi-lobes along the wires. In particular, 
the lowest boundary of the first MI lobes becomes pro- 
gressively flatter as a function of the hopping, inverting 
eventually its slope, leading to a re-entrant configuration 
(that remains in 2D). We have shown that such an effect 
leads to a non-trivial behavior of the MI plateaux in ex- 
periments with an axial harmonic confinement 0, Q . In 
particular, the MI plateaux may (for low hopping) be- 
come insensitive to the hopping, or even anti-intuitively 
grow for larger tunneling. Finally, we would like to stress 
that our results also apply to two-component Bose gases, 
predicting exciting phenomenology in on-going experi- 
ments in bosonic mixtures in lattices. 
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